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I. INTRODUCTION 



Multipole moments play an important role in Newtonian gravity and Maxwellian elec- 
trodynamics. Conceptually, there are two distinct notions of multipole moments — source 
multipoles which encode the distribution of mass (or charge-current), and field multipoles 
which arise as coefficients in the asymptotic expansions of fields. In Newtonian gravity, the 
first set is of direct interest to equations of motion of extended bodies while the second 
determines the gravitational potential outside sources. However, via field equations one can 
easily relate the two sets of multipoles. In the Maxwell theory, the rate of change of the 
source dipole moment is directly related to the energy fiux measured at infinity. Because 
of such useful properties, there has been considerable interest in extending these notions to 
general relativity. 

Results of the Maxwell theory were extended to the weak field regime of general rela- 
tivity — i.e., linearized gravity — quite some time ago. Already in 1916 Einstein obtained 
the celebrated formula relating the rate of change of the quadrupole moment of the source 
to the energy fiux at infinity [1]. In the fifties, Sachs and Bergmann extended the relation 
between the source multipoles and asymptotic fields [2]. However, in the framework of exact 
general relativity, progress has been slow. In the seventies, Geroch, Hansen and others [3] 
restricted themselves to the stationary context and introduced field multipoles by analyz- 
ing the asymptotic structure of suitable geometric fields constructed from the metric and 
equations they satisfy near spatial infinity. As in electrodynamics — and, in contrast to the 
situation in the Newtonian theory — they found that there are two sets of multipoles, the 
mass multipoles M(„) and the angular momentum multipoles J(„). In static situations, all 
the angular momentum multipoles J(„) vanish and the mass multipoles M(„) are constructed 
from the norm of the static Killing field which, like the Newtonian gravitational potential, 
satisfies a Laplace-type equation outside sources. In the (genuinely) stationary context, the 
J(„) are non-zero and are analogous to the magnetic multipoles in the Maxwell theory. In the 
Newtonian theory, since the field multipoles are defined as coefficients in the 1/r expansion 
of the gravitational potential, knowing all multipoles one can trivially reconstruct the po- 
tential outside sources. In general relativity, there are considerable coordinate ambiguities 
in performing asymptotic expansions of the metric. Therefore, Geroch and Hansen were 
led to define their multipoles using other techniques. Nonetheless, Bcig and Simon [4] have 
established that the knowledge of these multipoles does suffice to determine the space-time 
geometry near infinity. Their construction also shows that, if two stationary space-times 
have the same multipoles, they are isometric in a neighborhood of infinity. 

In the Geroch-Hansen framework, one works on the 3-dimensional manifold of orbits of 
the stationary Killing field and multipole moments M(„) and the J(„) arise as symmetric, 
trace-free tensors in the tangent space of the point A at spatial infinity of this manifold. 
Now, the vector space of n-th rank, symmetric traceless tensors on is naturally isomorphic 
with the vector space spanned by the linear combinations X]m=-n (^nmY^iO, (p) of spherical 
harmonics on the unit sphere in R^. Therefore, each M(„) and J(„) uniquely defines 
a set of numbers, M„„j and Jnm with m G {— n, — n + l,...,n — l,n}. Finally, let us 
consider stationary space-times which are also axi-symmetric (e.g., the Kerr space-time). 
Then Mn,m — MnSm,o and Jnm — Jn^mfl-i wheuce multipoles are completely characterized by 

two sets of numbers, M„, J„ with n = 0, 1, This fact will be useful to us because most 

of our analysis will be restricted to axi-symmetric isolated horizons. 

In stationary space-times, then, the situation with field multipoles is quite analogous to 
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that in the Newtonian theory. The status of source muhipoles, on the other hand, has not 
been so satisfactory. Dixon developed a framework to define source multipoles [5] but the 
program did not reach the degree of maturity enjoyed by the field multipoles. On the other 
hand, it is the source multipoles, rather than the field multipoles defined at infinity, that 
appear to be more directly relevant to the motion of extended bodies in general relativity. 
(For recent discussions, see [6].) In particular, it would be interesting to obtain useful 
generalizations of the quadrupole formula to fully relativistic objects such as neutron stars 
and black holes. Is there, for example, a relation between the flux of energy falhng across the 
horizon and the rate of change of the quadrupole moment of a black hole? To analyze such 
issues one needs a reliable definition of the source quadrupole moment, and more generally, 
of source multipoles. 

In this paper, we will focus on the problem of defining the analogs of source multipoles for 
black holes in equilibrium. The isolated horizon framework provides a suitable quasi-local 
arena for describing such black holes [7-11]. Thus, our task is to introduce an appropri- 
ate definition of multipoles which capture distortions of the horizon geometry and of the 
distribution of angular momentum currents on the horizon itself, and explore properties of 
these multipoles. Do these multipoles provide a diffeomorphism invariant characterization 
of the horizon geometry? If so, they may play a significant role in the analysis of equations 
of motion. Do they suffice to determine space-time geometry in the neighborhood of the 
horizon in stationary space-times? These are attractive possibilities. But since the hori- 
zon lies in a genuinely strong field region, a priori it is not obvious that a useful notion of 
multipoles can be introduced at all. Indeed, since the problem of defining and analyzing 
multipoles has turned out to be so difficult already for relativistic fluids, at first glance it 
may seem hopelessly difficult for black holes. However, a tremendous simplification occurs 
because black holes are purely geometric objects; one does not have to resolve the messy and 
difficult issues related to the details of matter sources. We will see that this simphfication 
makes it possible to carry out a detailed analysis and satisfactorily address the issues raised 
above. 

We will restrict our detailed analysis to axi- symmetric (or type II) isolated horizons, where 
the symmetry restriction applies only to the horizon geometry (and to the pull-back of the 
Maxwell field to the horizon) and not to the entire space-time. The material is organized as 
follows. In section II, we recall the relevant facts about isolated horizons and axi-symmetric 
geometric structures. In section III we define the multipoles M„ and J„ and show that these 
two sets of numbers provide a complete characterization of the isolated horizon geometry. 
For the specific case of Maxwell sources, we introduce another pair, P„ of electromagnetic 
multipoles and show that they suffice to determine the pull-back of the Maxwell tensor as 
well as that of its dual. Finally, the initial value formulation based on two intersecting null 
surfaces [12] implies that, if there is a stationary Killing field in the neighborhood of the 
horizon, multipoles also suffice to determine the near horizon geometry [13] . In section IV, 
we discuss potential extensions and applications of the framework ranging from numerical 
relativity to quantum gravity. 

Unless otherwise stated, in this paper all manifolds and fields will be assumed to be 
smooth. 
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II. PRELIMINARIES 



A. Isolated horizons 

In this sub-section we briefly recaU the relevant notions pertaining to isolated horizons 
[7, 9-11]. This discussion will also serve to fix our notation. 

Let us begin with the basic definitions [9]. A non- expanding horizon A is a null, 3- 
dimensional sub-manifold of the 4- dimensional space-time {A4,gab), with topology x M, 
such that: 

i) the expansion 6e of its null normal £ vanishes; and, 

ii) field equations hold on A with stress energy. Tab, satisfying the very weak requirement 
that —T'^b^'^ is a future-directed, causal vector. (Throughout, will be assumed to be future 
pointing.) 

Before discussing their consequences, let us note two facts about these assumptions: i) if 
the expansion vanishes for one null normal, it vanishes for all; and, ii) the condition on the 
stress energy is satisfied by all the standard matter fields provided they are minimally coupled 
to gravity. (Non-minimally coupled matter can be incorporated by a small modification of 
this condition. See, e.g., [14].) Since A is a null 3-surface, its intrinsic 'metric' qab has 
signature 0,+,+. The definition ensures that T^tf* oc and Ceqab = 0, where = denotes 
equality restricted to points of A. In particular, the area of any 2-sphere cross-section is 
constant on A. We will denote it by a a- The definition also implies that the space-time 
derivative operator V naturally induces a unique derivative operator D on A. Furthermore, 
DJ^ — uj^ for some globally defined 1-form cOa on the horizon. This 1-form will play an 
important role. 

The pair [q^b, D) is referred to as the intrinsic geometry of A. The notion that the black 
hole itself is in equilibrium is captured by requiring that this geometry is time independent: 
An isolated horizon (A, [£]) is a non-expanding horizon A equipped with an equivalence class 
of null normals where £"■ ~ £'^ if and only if £'"' — c£"' for a positive constant c, such that 
[Ci, D]=0 on A.^ 

Since A is a null surface, given any one null normal £"" in [£], we have £"'Da£^ = k^£^ for 
some K^. (Thus, — £"'U!a-) The requirement [Ci, D]^0 further implies that is constant 
on A. It is referred to as the surface gravity of A with respect to £"'. Note that if = c£" 
then K,£/ = CKf. Thus, the value of surface gravity refers to a specific null normal; it is 
not a property of (A, [£]). However, one can unambiguously say whether the given isolated 
horizon is extremal (i.e., has Kg, = 0) or non-extremal (i.e., has /t^ 7^ 0). 

The isolated horizon definition extracts from the notion of the Killing horizon just that 
'tiny' part which turns out to be essential for black hole mechanics [9, 11] and, more generally, 
to capture the notion that the horizon is in equilibrium, allowing for dynamical processes 
and radiation in the exterior region [7]. Indeed, Einstein's equations admit solutions with 
isolated horizons in which there is radiation arbitrarily close to the horizons [15]. Finally, 
note that the deflnition uses conditions which are local to A. Hence, unlike event horizons 
one does not require the knowledge of full space-time; the notion is not 'tclcological'. 

Of particular interest to our analysis is the case when the only matter present at the 
horizon is a Maxwell field. In this case, only the pull-hacks Bab and Eab to A, respectively of 



Given (A,ga5,-D), one can show that, generically, there is an unique equivalence class [f] of null normals 
such that (A, [£]) is an isolated horizon. However, our analysis will not be restricted to this case. 
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the Maxwell field Fab and its dual *Fab are needed in the isolated horizon analysis. We will 
refer to the quadruplet {qab, D, Bab, Eab) as the Einstein- Maxwell geometry of the isolated 
horizon. 

Remark: Note that the notion of an isolated horizon (A, [£]), can be formulated in- 
trinsically, using only those fields which define its geometry, without reference to the full 
space-time metric or curvature. Specifically: i) the condition is a null normal to A is 
captured in the property that qab has signature 0,-|-,-|- with qabi'' = 0; ii) the condition = 
can be replaced by Ctqab = 0; iii) the requirement on the stress-energy refers only to fields 
Bab and and iv) the condition [C^, D]^0 refers only to D (and not to the full space- 
time connection). This is why the quadruplet {qab, D, Bab, Eab) , defined intrinsically on A, 
was singled out to introduce the notion of Einstein-Maxwell horizon geometry. All equa- 
tions that are required in the derivation of the laws of mechanics of isolated horizons in 
Einstein-Maxwell theory [9, 11] as well as the main results of this paper refer just to these 
fields. 

Next, let us examine symmetry groups of isolated horizons. A symmetry of 
{Aji'^jqab, D, Bab, Eab) is a diffeomorphism on A which preserves qab, D, Bab and Eab and 
at most rescales i"' by a positive constant. It is clear that diffeomorphisms generated by 

are symmetries. So, the symmetry group Ga is at least 1-dimensional. The question 
is: Are there any other symmetries? At infinity, we generally have a universal symmetry 
group (such as the Poincare or the anti-de Sitter) because all metrics under consideration 
approach a fixed metric (Minkowskian or anti-de Sitter) there. In the case of the isolated 
horizons, generically we are in the strong field regime and space-time metrics do not ap- 
proach a universal metric. Therefore, the symmetry group is not universal. However, there 
are only three universality classes: 

i) Type I: the isolated horizon geometry is spherical; in this case, Ga is four dimensional; 

ii) Type II: the isolated horizon geometry is axi-symmetric; in this case, Ga is two dimen- 
sional; 

iii) Type III: the diffeomorphisms generated by ^" are the only symmetries; Ga is one di- 
mensional. 

Note that these symmetries refer only to the horizon geometry. The full space-time metric 
need not admit any isometrics even in a neighborhood of the horizon. Physically, type II 
horizons are the most interesting ones. They include the Kerr- Newman horizons as well as 
their generalizations incorporating distortions (due to exterior matter or other black holes) 
and hair. Our main results refer to the type II case. 

Finally, one can use the field equations to isolate the freely specifiable data which de- 
termines the Einstein-Maxwell geometry of an isolated horizon [10, 16]. The analysis is 
naturally divided in to two cases: Ke ^ and K£ = 0. Fix any 2-sphere cross-section A of A 
and denote by q^-b the natural projection operator on A. Then, on any type II horizon, we 
have the following.^ 

• In the non-extremal case, the free data consists of the projections qab, oJa, Bab, Eab to A 
(using q"-b) of qab,'^a, Bab, Eab on A. (Recall that uja is defined by Dad-^ = (jJa^^-) That 
is, given the free data on A, using the projections of the field equations on A, one can 



^ For convenience of the reader, proofs of these assertions are sketched in Appendix A. Throughout this 
paper, the term 'free data' refers to freely specifiable fields. In order to rcconstriict the horizon geometry, 
in addition to specifying this data one also needs to fix some constants, such as the horizon area oa- 
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reconstruct the full Einstein-Maxwell geometry {qab, D, Bab, Eab) on A. 

• In the eoctremal case, remarkably, the fields qab, ijJa, Bab and Eab turn out to be all 
universal; they are the same as those in the extremal Kcrr-Newman geometry [16]! The 
free data consists of {Sab), a symmetric, second rank tensor on A, determined by D via 
Sab = q.aJlb'^Dcnd where Ua is the covariant normal to A within A, satisfying £"na = — 1. 
(Thus, Sab is the analog of the 'extrinsic curvature of A within A.) Given a Kerr- 
Newman quadruplet qab, oJa, Bab and Eab the only freedom in the Einstein-Maxwell 
horizon geometry lies in the part of D which determines the 'extrinsic curvature' Sab 
of A. 

These results will play an important role in section 3. Specifically, we will construct multi- 
poles from the free data and show that they provide a diffeomorphism invariant characteri- 
zation of the free data and hence, via field equations, of the horizon geometry.^ 



B. Axi-symmetric structures 

In this sub-section we will recall a few facts about axi-symmetric geometries on S^. This 
discussion will also make our assumptions explicit. 

Let S be a manifold with the topology of a 2-sphere, equipped with a metric qab- We will 
denote by e<jb the alternating tensor on S compatible with qab, by a the area of S and by R 
its radius, defined by a = AttR'^. We will say that (S, qab) is axi-symmetric if it admits a 
Killing field 0" with closed orbits which vanishes exactly at two points of S. The two points 
will be referred to as poles. We will now show that such metric manifolds carry invariantly 
defined coordinates and discuss properties of the metric coefficients in these coordinates. 

Since C^f^eab — 0, there exists a unique (globally defined) function C on S such that 

£'aC = ;^4a/ and j^Ce = 0. (2.1) 

It is clear that = and DaC vanishes only at the poles. Hence, C is a monotonic function 
on the 1-dimensional manifold § of orbits of 0". Now S is a closed interval and the end points 
correspond to the two poles. Hence C monotonically increases from one pole to another. We 
will say that C assumes its minimum value at the south pole and the maximum at the north 
pole. 

Next, let us introduce a vector field on S' = (S — poles) via: 

qabCy'^O and 1 (2-2) 

Then it follows that 

C = ^r''DbC, (2.3) 

where /x^ = qab4>"'<P'^ is the squared norm of 0". Hence integral curves of go from the 
south pole to the north (and diverges as one approaches the poles). Using C", we can 



^ This general philosophy is rather similar to that used by Janis and Newman to introduce a notion of 
multipoles using 'free data' on null surfaces [17]. However, since our analysis is restricted to isolated 
horizons rather than general null surfaces, our results are free of coordinate ambiguities. 
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define a preferred affine parameter of 0" as follows: Fix any one integral curve / of 
in §' and set = on /. Thus, on /, we have £^0 = 0. Now, since 0" is a Killing field 
and is constructed uniquely from {qab, 0"), it follows that £,/,C" = 0. Hence, we conclude: 
C^{Cq4>) — 0. Since £^0 = on 7, it now follows that — const on every orbit of C"- This 
now implies that the affine parameter of 0° has the same range on every orbit of 0".^ 
Without loss of generality we will assume that 0" is such that G [0, 2ti) on S. 

Thus, starting from geometry, we have constructed two coordinates on § such that 
0" = {d/d(j)Y and = {d/dQ"' are orthogonal. Eqs (2.1) and (2.3) now imply that the 
metric has the form: 

qab^R\r'D,CD,C + fDa(t>D,(t>) and g-^^ = -1 (/(«('' + /"W) (2-4) 

where / = \j? jB?. The fact that the area of S is ^'kB? now imphes that the range of C is 
necessarily [—1, 1]. Conversely, given any 'coordinates' C', 4>' in [—1, 1] x '^I'^.'k'L in which the 
metric can be expressed in (the primed version of) the form (2.4), it is easy to show that 
C' = C 3^nd 0' is at most a rigid shift of 0. (For an analogous construction in Newman-Penrose 
notation, see [18].) 

Functions C, serve as 'coordinates' on S modulo usual caveats: they are ill-defined at 
the poles and has a 27r discontinuity on the integral curve / of C"- We have to ensure 
that the metric qah is smooth in spite of these coordinate problems. The discontinuity at / 
causes no problems. However, poles do require a careful treatment because the norm of 
0" — and hence / — vanishes there. Smoothness of qah at the poles (i.e. absence of conical 
singularities) imposes a non-trivial condition on /: 

^hm /'(C) = T2 (2.5) 

where 'prime' denotes derivative with respect to C,. On a metric 2-sphere, we have / = 1 — 
and we can bring the metric to the standard form simply by setting C = cos Q. In the general 

case, / has the same values and first derivatives at the poles as on a metric 2-sphere. Using 
I'Hopital's rule, one can show that this fact suffices to ensure that the metric (2.4) is smooth 
at the poles. 

Finally, we note a property of these axi-symmetric metrics which will be useful in section 
HI B. A simple calculation shows that the scalar curvature of qah is given by: 

nC^A)^-^J'\Q. (2.6) 

By integrating it twice with respect to C, and using the boundary conditions (/|^=_i) ~ 
and (/'|^=_i) = 2, one can reconstruct the function / from the scalar curvature: 



^ ~ ■ +2(C + 1) (2.7) 



dCi /'dC2^(C2) 
1 J -I 



Thus, thanks to the preferred coordinates adrnitted by an axi-symmetric geometry on S, 
given the area a of S and the scalar curvature TZ, the metric qah is completely determined. 



^ To make this argument precise, one should work on the universal covering of the orbits of (jf" on S' but 
the additional steps are straightforward. 
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Remark: Coordinates {C,<P), determined by the axi-symmetry of Qab, also enable us to 
define a canonical round, 2-sphere metric on §: 

= RVo'DaC D,C + foDaCP D,ct>) (2.8) 

where fo = 1 — C^. Note that q^^ has the same area element as qab- This round metric 
captures the extra structure made available by axi-symmetry in a coordinate invariant way. 
The availability of enables one to perform a natural spherical harmonic decomposition 
on §. This fact will play a key role throughout section III. 

III. MULTIPOLES OF TYPE II ISOLATED HORIZONS 

This section is divided into four parts. In the first three, we restrict ourselves to non- 
extremal, type II isolated horizons with no matter fields on them. We begin in the first part 
by defining a set of multipoles, /„, L„, starting from the horizon geometry. We then show 
that one can reconstruct the horizon geometry starting from these two sets of numbers. In 
the second, we show that if two isolated horizons have the same multipoles /„, L„, their 
geometries are related by a diffeomorphism. Thus, these two sets of numbers provide a 
convenient diffeomorphism invariant characterization of the horizon geometry. Therefore, 
we refer to /„, L„ as geometric m,ultipoles. In the third part, we rescale these moments by 
appropriate dimensionful factors to obtain mass and angular momentum multipoles M„, J„. 
Finally, in the fourth part we first discuss electromagnetic multipoles and then summarize 
the situation on extremal isolated horizons. 

A. Geometric Multipoles 

Let (A, [£]) be a non-extremal, type II isolated horizon with an axial Killing field 0". 
In this sub-section, we will ignore matter fields on A and concentrate just on the horizon 
geometry defined by {qab,D)- Fix a cross-section A of A. Then, as summarized in section 
II A, the free data that determine the horizon geometry consists of the pair {qab,i^a) where 
(Jab is the intrinsic metric on A and Ua is the projection on A of the 1-form Ua on A (defined 
by Dai'^ = ujaf-^)-^ However, there is some gauge freedom associated with our choice of 
the cross-section A [10]. We will first spell it out and then define multipoles using gauge 
invariant fields. 

For simplicity of presentation, let us fix a null normal in Then, A can be regarded 
as a leaf of a foliation u = const such that £°'DaU = 1. For notational simplicity, we will set 
Ua = —DaU so that ria is the covariant normal to the foliation satisfying I'^ria = — 1. The 
projection operator q^ on the leaves of this foliation is given by qa = 5^ + riad'^- Hence, 
Q.ab — Qab ^ud uja — i^a-\- i^i^a as tcusor fields on A. Since Cf^qah = and Cf^uja = on any 
isolated horizon, and since Cgria^Q from the definition of rta-, it follows that {qabi^a) on any 
one leaf is mapped to that on any other leaf under the natural diffeomorphism (generated 
by relating them. Let us now consider a cross-section A' which does not belong to this 
foliation. Let u' — const denote the corresponding foliation. Set F — u — u' . Then, regarded 
as tensor fields on A the two sets of free data are related by 

Qab = Qab and Oj'^ = UJa + lit DaF (3.1) 
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ThuSj^ under the natural diffeomorphism (defined by the integral curves of between A 
and A', qab is mapped to but Ua is not mapped to a)^; the difference is a gradient of a 
function. This is the gauge freedom in the free data. 

It is therefore natural to consider, in place of its curl. Prom the isolated horizon 
framework, it is known that for any null tetrad £",n",m",m" such that e the Weyl 
components ^I/q and vanish on A whence ^2 is gauge invariant [9, 10], and the curl of uja 
is given just by Im^2: 

^[aCJt] = Im^2 Cafe (3.2) 

where eab is the natural area element on A (which satisfies = and C(,eab = ^)- Thus, 
the gauge invariant content of Co a is coded in Im^2- 

The second piece of free data is the metric qah on A. In section II B we showed that 
using an invariant coordinate system (C,0), one can completely determine qah in terms of a 
number, the area a and a function, its scalar curvature TZ. If, as assumed in this sub-section, 
the cosmological constant is zero and tliere are no matter fields on A, then TZ = — 4Re\l'2 
[8]. Hence the gauge invariant part of the free data that determines the horizon geometry 
is neatly coded in the Weyl component ^2- 

It is therefore natural to define multipoles using a complex function $a on A: 

$^ := i:^-ilm^2. (3.3) 
4 

(Thus, in absence of matter on A, $a = —^2 while in presence of matter it is given by 
$A = — ^2 + (l/4)-Ra69''''— (1/12)-R, where Rah is the Ricci tensor and R the scalar curvature of 
the 4-metric at the horizon.) Since all fields are axi-symmetric, using the natural coordinate 
C on A, we are led to define multipoles as: 

7„ + iL„:= jf $Ai;°(C)d'V^ (3.4) 

or, 

4:=^£^F„°(C)d2y and K := - ^^Im^^Y^iQ d'V . 
Here are the m — spherical harmonics, subject to the standard normalization: 

^jXd'V^Rl5r,,m: (3-5) 

where i?A is the horizon radius defined through its area a a via a a = 4:7iR\. Thus, given 
any horizon geometry, we can define a set of two numbers, and L„. Recall that C is 
determined entirely by the metric qab and the rotational Killing field 0". Therefore, it is 
immediate that if (A, qab, D) and (A', g^^, D') are related by a diffeomorphism, we have 
In = I'n and Ln = L[^\ the two sets of numbers are diffeomorphism invariant. If the isolated 
horizon were of type I, qab would be spherically symmetric and Im^'2 would vanish [9]. 
Then, the only non-zero multipole would be /q which, by the Gauss-Bonnet theorem, has a 
universal value (see below). Given a generic type II horizon, as we saw in section II B, the 
axi-symmetric structure provides a canonical 2-sphere metric (which, in the type I case, 
coincides with the physical metric). The physical geometry has distortion and rotation built 
in it. The round metric q^h serves as an invariantly defined 'background' against which one 
can measure distortions and rotations. Multipoles /„, L„ provide a diffeomorphism invariant 
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characterization of these. More precisely, they encode the difference between the physical 
horizon geometry {qab,D) and the fiducial, type I geometry determined by (g°5,u;a = 0). 

Finally, note that In and L„ can not be specified entirely freely but are subject to certain 
algebraic constraints. The first comes from the Gauss-Bonnet theorem which, in the axi- 
symmetric case, follows from the boundary condition (2.5) on /' and the expression (2.6) of 
the scalar curvature in terms of /: 

/o = ^ ^ro°(C)d'V = x/^. (3.6) 

The second comes directly from the relation (3.2) between Im\E'2 and curl oi uja'- 

Lo = - ^= / Im*2 d^y = 
V47r 7 a 

The third constraint comes again from (2.5) and (2.6): 

8V^ /a 

We will show in section IIIC that (3.7) implies that, as one would physically expect, the 
'angular momentum monopole' necessarily vanishes and (3.8) implies that the mass dipole 
vanishes, i.e., that our framework has automatically placed us in the 'center of mass frame 
of the horizon'. Next, because $a is smooth, these moments have a certain fall-off. Let us 
assume that $a is (i.e., the space-time metric is C*"*"^). Then as n tends to infinity, /„ 
and Ln must fall off in such a way that 

oo fc+1 

5] ^|n'"(4-iL„)|2<oo (3.9) 

n=0 m=0 

Finally, there is a constraint arising from the fact that / is non-negative and vanishes only 
at the poles. (This property of / is essential for regularity of the metric.) Using (2.6) and 
the definition (3.4) of /„, one can express / in terms of /„. Unfortunately, the resulting 
restriction on multipoles is quite complicated: 

oo 

^ Vvr(2n + 1) 2n + 3 

2(2n+l) IQU >o (3.10) 

(2n + 3)(2n-l) "^^^ 2n - 1 " ^^^^J " - ^ > 

and vanishes only if = ±1, where P„ are the Legendre polynomials. 

Remark: We conclude by noting some simplifications that occur in presence of additional 
symmetries. Certain space- time metrics, such as the Kerr solutions have a discrete (spatial) 
reflection symmetry, (C, 0) i— > (— 4>'^^)-i under which ^2 ^ ^2- Therefore, in the isolated 
horizon framework it is interesting to consider the case in which $a ^ $a ^'^d 1-^ Eab 
under the discrete diffeomorphism ^ —( on A. Then, since are even/odd under 
reflections if n is even/odd, it follows that /n = for all odd n and L„ = for all even n. 
Next, consider the case in which the isolated horizon is a Killing horizon of a static Killing 
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(3.7) 



(3.8) 



field. Then, one can show that Im^2 = [11]. Hence L„ = for all n but in general /„ can 
be arbitrary, capturing possible distortions in the horizon geometry. Finally, consider type I 
isolated horizons on which the horizon geometry is spherically symmetric. Then, $a and eab 
are spherically symmetric. It is then obvious from properties of spherical harmonics Y^{() 
that for all n > 0, wc have In — — 0. Since Lq always vanishes, in this case the only 
non-trivial multipole is Iq — \pK. 



B. Reconstruction of the horizon geometry and the uniqueness issue 

1. Construction 

Let us first show that the knowledge of the area and multipole moments suffices to 

reconstruct an isolated horizon geometry. 

Suppose we are given only the radius and the set {In, L^} of geometric multipoles of 
a , non- extremal, type II isolated horizon. Then, one can explicitly construct a , non- 
extremal, type II isolated horizon geometry (A, qah, D) such that area is given by oa = 47ri?^, 
and geometric multipoles are given by In and Ln- 

Let S be a smooth 2-manifold, topologically S^. Let (C,0) be a spherical coordinate 
system on §, with ^ G [—1, 1] and (p ^ [0, 27r). Define 

A °° 

^ ■= wT^^nY^ (3.11) 



Im^-; 



^ oo 

-wT^^nY^iO- (3.12) 



A n=0 



Since {/„, L„} are multipoles of a isolated horizon, we are guaranteed that the fields TZ 
and Im^'2 so defined are C^. Next, define a function / on § via: 



/(C) := 

and set 



^^1^(1 y"%c2^(C2) 



2(C + 1) (3.13) 



Qab Rl ^aCAC + /(C) ^a^A^j ■ (3.14) 

Since arc the multipoles of an isolated horizon, (3.10) is satisfied whence it follows that / 
is non-negative on § and vanishes only at the 'poles', ^ = ±1. Hence Qab is a smooth metric 
except possibly at the poles. Next, /q = \/ti and /i = 0. Interestingly, these conditions 
imply /'(±1) = ^2, and therefore ensure that the metric Qab can be smoothly extended to 
points ( = ±1 oi E>. Next, it is straightforward to verify that the area of S with respect to 
this metric is given by AnR^ and its scalar curvature is given by IZ. Let us now turn to the 
multipoles L„. Since these multipoles also come from a type II isolated horizon geometry, 
we have Lq = 0. Eq (3.12) now imphes f Im^2 d'^V — 0, whence there is a globally defined 
1-form Ua on § such that 

DiaUb] = Ini*2 ^ab (3.15) 
where eab is the alternating tensor on (S, Qab) ■ 
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Now consider a 3-manifold A = § x M and equip it with a vector field along the 
'M-direction'. Let vr be an embedding of S onto a 2-sphere cross-section A of A. Then, the 
diffeomorphism tt equips A with fields qab, i^a- This is the free data required to construct a 
non-extremal isolated horizon geometry {qab,D) on A (see [10] or Appendix A). Since the 
data are axi-symmetric, so is the horizon geometry. Finally, from our definitions of multipole 
moments, it follows immediately that the multipole moments of this isolated horizon are 
given by {/„,!/„}. 

Remark: Our construction in fact suffices to show the following stronger existence result: 
Given a positive number i?A and a set of real numbers /„, L„ subject to the constraints 
Eqs (3.6)-(3.10), there is a non-extremal isolated horizon geometry {A,qab,D) with area 
OA = 47ri?^ and geometric multipoles /„, L„. We chose not to refer to the conditions (3.6)- 
(3.10) but assumed instead that the /„, L„ arise from a type II isolated horizon for simplicity. 
Although the existence result is stronger, it is difficult to verify if a given set (/„, L„) satisfies 
Eqs (3.9) and especially (3.10). The overall situation is similar at spatial infinity. There, 
the constraints that must be imposed on field multipoles to ensure that the formal power 
series converges to a non-degenerate metric are awkward to state and investigations have 
focussed on reconstruction [3] . 

2. Uniqueness 

The question now is whether the geometry we constructed in section III B 1 is diffeomor- 
phic to the one from which the multipoles were first constructed. 

Let (A, [i], Qab, D) and (A', [£'], g^;,, D') he two type II non-extremal isolated horizons with- 
out matter fields, with same area and same geometric multipoles In,Ln- Then there is a 
diffeomorphism from A' to A which maps {[£'], q'^^, D') to {[£],qab, D)- 

Choose null normals i G [£] and £' G [£'] such that K,i = kii = k and consider foliations 
u = const and u' = const of A and A' which are tangential to 0" and 0'" respectively such 
that £°'DaU — 1 and £"^Dau' — 1. The intrinsic metrics on these leaves are axi-symmetric. 
Hence, as in section II B, we can find invariant coordinates (C, 0) and (C', (p') on the two sets 
of cross-sections. Let vr be the diffeomorphism from A to A' defined by {u, (, (j)) ^ {u', (f)'). 
Since the multipoles /„ are the same, it follows from (3.11) that vr maps TZ' to TZ and hence 

^'^ 'lab- Similarly, since the multipoles L„ are the same, it follows from (3.12) that n maps 
Im\E'2 on A' to Im^2 on A. But this only implies that da;' is mapped to duj, i.e. 

TT* Oj'^^0a + nDah (3.16) 

for some function h. Thus, in general, tt does not map the free data on A' to that on A. 

However, as the discussion in the beginning of sub-section III A suggests, this can be easily 
remedied by changing the foliation on A. Define u = u — h on A. Then, the projection 
LJa of uJa ou the ucw foliatiou is given by LJa = + nDah. Hence the diffeomorphism 
7f : {u, C, 0) I— ^ {u', C', 0') from A to A' does map the primed free data to the barred free 
data: 

Tt* u'^ = Qa and 7f* q'^ = qab = qab (3.17) 

This isomorphism between the two sets of free data naturally extends to an isomorphism 
between the two horizon geometries (see [10] or Appendix A). 

Thus, together with area, multipoles /„, L„ provide a convenient, diffeomorphism invari- 
ant characterization of type II horizon geometries. 
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C. Mass and angular Momentum multipoles 

As is obvious from their definition, L„ are all dimensionless. Therefore, it is difficult 
to attribute a direct physical interpretation to them. In this sub-section we will argue that 
they can be rescaled in a natural fashion to obtain quantities which can be interpreted as 
mass and angular momentum multipoles M„ and J„. 

In the isolated horizon framework, the area oa is defined geometrically. One then defines 
the horizon angular momentum Ja as a surface term in the expression of the generator of 
rotations, evaluated on a 2-sphere cross-section of the horizon. Ja is unambiguous because 
type II horizons come with an axial symmetry [11]. The horizon mass Ma is also defined 
using Hamiltonian methods as the generator of a preferred time translation. However, the 
preferred time translation varies from space-time to space-time. If Ja = 0, it points along 
if not, it is a suitable linear combination of and 0", which can be fixed only after 
one knows the value of Ja [11]. Thus, Ja is defined first before one can fix Ma. In the 
same spirit, we will first define the angular momentum multipoles J„ and then the mass 
multipoles M„. 

We begin by recalling a general fact about angular momentum. Fix a space-time {A4, gab) 
and a space-like 2-sphere S in it. Let (p be any vector field tangential to S. Then, by regard- 
ing S as the inner boundary of a partial Cauchy surface M, one can use the Hamiltonian 
framework to define a 'conserved' quantity Jg 



where Kab is the extrinsic curvature of M. In a general space-time, this quantity is inde- 
pendent of M if and only if (p is divergence free with respect to the natural area element 
of S. Thus, for each divergence-free (/? on S, Jg depends only on S and can be interpreted 
as the (^"-component of a 'generalized angular momentum' associated with S. If S happens 
to be a cross-section of A, as one would expect, one can recast this expression in terms of 
the fields defined by the isolated horizon geometry, making no reference at all to the partial 
Cauchy surface M [11]: 



Here / is a 'potential' for (p"- on A — given by = e°'^Dbf — which exists because C^eab = 0. 
By the isolated horizon boundary conditions it follows that if is the restriction to of a 
vector field on A satisfying Ce, pf" = 0, then Jl" is independent of the 2-sphere cross-section 
S used in (3.19). 

Thus, on any isolated horizon there is a well-defined notion of a 'generalized angular 
momentum' J|', associated with any divergence free vector field 99" satisfying Liipf - = 0. 
Im^'2 plays the role of the 'angular momentum aspect'. Hence, it is natural to construct the 
angular momentum multipoles J„ by rescaling the L„ with appropriate dimensionful factors. 
This strategy is supported also by other considerations. First, since Im^2 transforms as 
a pseudo-scalar under spatial refiections, we will automatically satisfy the criterion that 
the angular momentum multipoles should transform as pseudo tensors. Second, all angular 
momentum multipoles would vanish if and only if Im^'2 = and this is precisely the condition 
defining non- rotating isolated horizons [9, 11]. Thus, the strategy has an overall coherence. 




(3.18) 




(3.19) 
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To obtain the precise expression, let us first recall the situation in magnetostatics in flat 
space-time. If the current distribution is axi-symmetric, the nth magnetic moment M„ 
is given by: 

J r''Pn{cos9)V ■ {xxj)d^x, (3.20) 

where P„ are the Legendre polynomials. If the current distribution is concentrated on the 
sphere S defined hy r — R, the expression simplifies to: 



M„ = -i^'^+i i {e^'DbPnicosO)) jad'V, (3.21) 



Is 

where iab is the alternating tensor on the r = R 2-sphere and jb is the projection of jb on 
this 2-sphere. Note that this expression refers just to the axi-symmetric structure on the 
2-sphere S and not to the flat space in which it is embedded. Comparison of (3.19) with 
(3.21) suggests that we can think of the horizon A as being endowed with a surface 'current 
density' 

(jA)a=^^a (3.22) 

and define the angular momentum (or 'current') moments as: 



Jn = j( (?"£l6P„(C))'i.<<V 



(3.23) 



2n + l AnG 

(3.24) 

Let us now turn to the mass multipoles, M^- When all J„ vanish, we should be left 
just with M„. These are then to be obtained by rescaling the multipole moments /„ by 
appropriate dimensionful factors. In electrostatics, when the charge density is axi-symmetric, 
the electric multipoles are defined by 

En= [ r''Pnicos9)pd^x. (3.25) 



i„ = i?" ^ p„(cos 9)p (3.26) 



When the charge is concentrated on the sphere S defined hy r — R, the expression simphfies 
to: 

Is 

where p is the surface charge density. Again, the final expression refers only to the axi- 
symmetric structure on the 2-sphere S and not to the fiat space in which it is embedded. 
Hence we can take it over to type II horizons. What we need is a notion of a 'surface mass 
density'. Now, Hamiltonian methods have provided a precise definition of mass Ma of type 
II isolated horizons in the Einstein- Maxwell theory [9, 11]. The structure of geometric mul- 
tipoles /„ now suggests that we regard Ma as being 'spread out' on the horizon, the 'surface 
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density' pa being uniformly distributed in the spherical case but unevenly distributed if the 
horizon is distorted. It is then natural to set^ 



= J- MiJl = — ^- MARe*2 ■ 



(3.27) 



This heuristic picture motivates the following definitions: 



V 2n + 1 27r 




i;°(C)Re*2ciV 



47r M^R\ ^ 



n • 



(3.28) 



2n + 1 27r 



Here Ma is the isolated horizon mass, which is determined by the horizon radius R/\ and 
angular momentum Ji via: 



Since our definitions are based on analogies with source multipoles in the Maxwell theory, 
an important question is whether they have the physical properties we expect in general 
relativity. We have the following: 

• As discussed in section III A, the geometrical multipoles Lq and /i vanish. Hence 
it follows that the angular momentum monopole moment Jq vanishes as one would 
expect on physical grounds, and the mass dipole moment Mi vanishes implying that 
we are in the center of mass frame. 

• By construction, the mass monopole Mq agrees with the horizon mass Ma and, by 
inspection, the angular momentum dipole moment Ji equals the horizon angular mo- 
mentum Ja, calculated through Hamiltonian analysis [11]. 

• Restrictions on the multipoles imposed by symmetries of the horizon geometry follow 
immediately from the remark at the end of section HI A: 

i) if the horizon geometry is such that i— > ^'2 under the reflection C, ^ — C, then 
M„ = if n is odd and J„ = if n is even. This is in particular the case for the Kerr 
isolated horizon. 

ii) If Im\l/2 = 0, then all angular momentum multipoles vanish. This is in particular 
the case if A is a Killing horizon in a static space-time. 

iii) If the horizon geometry is spherically symmetric, M„ = for all n > and J„ = 
for all n. 

• Our multipoles M„, J„ are constructed just from the knowledge of the horizon geom- 
etry; knowledge of the space-time metric in the exterior region is not required. In 
particular, there may well be matter sources outside the horizon, responsible for its 
distortion and the exterior geometry need not even be stationary or asymptotically 



^ Note incidentally that the smaller the principal radii of curvature of the intrinsic geometry, the higher is 
Pa- Thus, the situation has a qualitative similarity with the way charge is distributed on the surface of a 




(3.29) 



conductor. 
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flat. Even when the exterior is stationary and asymptoticafly flat, there is no a priori 
reason to expect that these 'source multipoles' would agree with the 'fleld muhipoles' 
deflned at inflnity — unless symmetry principles are involved — because the gravita- 
tional fleld outside the horizon would also act as a source, contributing to the 'total' 
moments at inflnity. 

It is interesting to analyze the situation in Kerr space-times. Because of symmetries, 
the horizon mass monopole and angular momentum dipole agree exactly with the 
corresponding fleld moments at inflnity. What about the mass quadrupolc M2 or the 
angular momentum octupole? As flgure 1 shows, while the fleld quadrupolc is not 
exactly equal to the horizon quadrupole, the difference is insigniflcant for small values 
of the Kerr parameter a = J/M . Furthermore, even in the extreme limit a = M, there 
is only a factor of 1.4 between the two mass quadrupoles and 1.14 between the angular 
momentum octupoles. In the large a limit, the effects due to dragging of inertial 
frames is so large that post-Newtonian and post-Minkowskian approximations fail; in 
this regime, there are no reliable calculations relating the source and fleld moments. 
So a priori it would not have been surprising if there were a factor of a million! 




FIG. 1: Plots of the horizon multipoles (solid lines) and Hansen's field multipoles (dashed lines) 
for Kerr space-times as functions of a/M. The first figure shows the behavior of |M2|/M^ where 
M2 is the mass quadrupolc moment and the second shows the behavior of \J^\/M^ where J3 is 
the angular momentum octupole. Absolute values are used because of the difference in Hansen's 
and our sign conventions. (His mass monopole is —M while ours is M.) Quantities plotted are 
dimensionless in the G = 1 units. 

• As noted in the Introduction, in vacuum, stationary space-times Hansen's fleld mo- 
ments suffice to determine the geometry in the neighborhood of inflnity. Is there an 
analogous result for the horizon multipoles? The answer turns out to be affirmative. 
Fix a cross section 5" of A and consider the future directed, inward pointing null vector 
fleld which is orthogonal to S, normalized so that i'^Ua = — 1. The null geodesies 
originating on S with tangent vector —n"- generate a null surface A/". In the source-free 
Einstein theory, there is a well-deflned initial value problem based on the double null 
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surfaces, A and H [12]. The freely specifiable data consists of the pair {qahi'^a) on A 
and the Newman- Penrose component of the Weyl tensor on M [7]. 

Now suppose that the space-time is analytic near A and admits a stationary Killing 
field t"' in the neighborhood of A which is time-like in the exterior region and be- 
comes null on A. Then, as one might imagine, the initial value problem becomes 
highly constrained: ^4 on H is determined by the horizon geometry [13]. Hence, the 
horizon multipoles suffice to determine the solution to Einstein's equation in a past 
neighborhood of A U N . Thus, qualitatively the result is the same as with the field 
multipoles.^ 

To summarize, the mass and angular momentum multipoles M„, J„ have physically ex- 
pected properties. This in turn strengthens the heuristic picture we used to fix the dimen- 
sionful rescalings of In,Ln. We first considered stationary, axi-symmetric charge and current 
distributions with support on a 2-sphere in the Maxwell theory and expressed the electric 
and magnetic multipoles using only the axi-symmetric structure on the 2-spheres without 
reference to Minkowski space-time. We then noted that these structures are available also on 
type 11 horizons.^ Additional structures made available on these horizons by geometric and 
Hamiltonian methods then led us to our definitions of M„, J„. The physical picture is that 
observers in the exterior region (between the horizon and infinity) can regard the horizon 
multipoles as arising from an effective (but fictitious) mass density Pa = — (l/27r)MARe\I'2 
and a current density (jA)a = (l/87rG)a;a on A. This picture may be useful in physical 
applications. 

D. Extensions 

So far, we have restricted ourselves to non-extremal, type II isolated horizons which have 
no matter fields on them. In this sub-section, we will extend our results in various directions. 

1. Maxwell fields 

Let us restrict ourselves to the non-extremal case but allow Maxwell fields on A. Then the 
Einstein-Maxwell horizon geometry consists of the quadruplet {qab, D, Bab, Eab)- The pres- 
ence of matter fields on A causes a few minor modifications in the discussion of gravitational 
multipoles. We will first discuss these and then turn to the electromagnetic multipoles. 



^ However, while the domain in which the solution is determined by field multipoles is known to be large 
[4] , the available results on the double- null initial value problem only ensure the existence of the solution 
in a small, past neighborhood of AuAf. But it is quite possible that the double-null initial value results 

can be significantly strengthened. 

Note that the spherical harmonics Y^{() used in the definition of geometric multipoles /„, Ln refer to the 
unique round metric q°f^ defined by coordinates (C,'/') (see (2.8)). They are eigenstates of the Laplacian 
defined by and not of the physical metric qab- Therefore, the extension of the Maxwell formulas is 
natural. As noted in section II B, q°f^ is uniquely determined by the axi-symmetry of qab and the two 
metrics have the same area element. Therefore, pA and (jA)a can be interpreted as the mass and current 
densities in terms of either metric. 
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In the gravitational case, the definition of the basic complex field $a on A continues to 
be the same in terms of TZ and Im^2 but changes if we use Re^2 in place of TZ: 

$A := - ilra^2 = -^2 + ^Rabq"''. (3.30) 

where the term involving the space-time Ricci tensor can be expressed in terms of Bab and 
Eab as: 

Rahr' = G {BabB'^' + EabE'^') . (3.31) 
The definition of mass Ma changes from (3.29) to 

2GR-A • ^^-^^^ 

The subsequent definition of the geometric and physical multipole moments is the same as 
in section III A. The reconstruction of the free data from multipoles is also unaffected but 
there is an additional term involving Bab and Eab in the reconstruction of D from the free 
data (see Appendix A). Finally, there is a minor modification in the list of properties of 
Mn, Jn listed in section IIIC. This arises because, in presence of a Maxwell field on the 
horizon, the canonical angular momentum Ja obtained by Hamiltonian methods contains 
two terms, a gravitational one and an electromagnetic one [11]. The angular momentum 
dipole moment Ji yields just 'the gravitational part' of Ja. There does not seem to be a 
natural generalization of the definition of angular momentum multipoles such that Ji would 
agree with the full Ja. 

Let us now turn to the electromagnetic fields. As noted in section II A, the electromag- 
netic free data consists of the projections Bab and Eab of Bab and Eab on a cross-section A. 
Therefore, it is natural to define the electromagnetic counterpart of $a: 

^T----l^^'[Eab + iBab] (3.33) 
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and define multipoles via: 



Q„ = 6 Re($™)i;°(C)dV (3.34) 



y/4:7r{2n + 1) A 

i?" r 

v/47r(2n + 1) Ja 



P„ = i lm{<^f ') F„°(C) d'V . (3.35) 

(3.36) 



Clearly, Qq and Pq are the electric and magnetic charges of the horizon. Thus, heuristically 
Re($™)/47r and Im($™)/47r may be thought of as 'surface charge densities' on the hori- 
zon and charge multipoles capture the non-uniformity in the distributions of electric and 
magnetic charge densities. 

Remark: If non-electromagnetic sources are also present, we can still define the gravita- 
tional (and electromagnetic) multipoles as above but the multipoles for other sources have 
to be defined case by case. Gravitational multipoles again determine the 'free data' for the 
horizon geometry {qab,D). However, to reconstruct the horizon geometry from this data, 
one needs to know those matter fields which determine the components Qa'^ Qb'^ Red of the 
Ricci tensor, evaluated at the horizon. 
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2. Extremal horizons 



Let us now summarize the situation with extremal horizons. In this case, important 
modifications are needed in the gravitational sector. As noted in section II A, surprisingly, 
{Qab,^a, Bab, Eab) ^^e uow Universal, determined by the Kerr-Newman family [16]. How- 
ever, they no longer suffice to determine the derivative operator D on the horizon because 
the second rank, symmetric tensor field Sab = Qa'^Qb'^ Dji^ is now free [10] (see Appendix 
A). Therefore, multipole moments characterizing the free data in the geometry have to be 
constructed from components of Sab- Again, as with uja, there is a gauge freedom in Sab 
associated with the choice of the cross-section A used to evaluate it. This can be eliminated 
by fixing the trace of Sab [13]. Then the gauge invariant information is coded in the two 
components of the trace-free part of Sab- 

Of course, we can still use $a to introduce the geometric multipoles /„, L„ and elec- 
tromagnetic multipoles Qn, P„. But they are universal; the same as in the extremal Kerr- 
Newman case. However, now there are two additional sets of geometric multipoles obtained 
by integrating the two independent components of Sab against Y^{Q. They capture the free 
data which can vary from one extremal isolated horizon to another. To our knowledge, they 
do not have a simple, direct physical interpretation. 

3. Type III horizons 

The notion of the horizon geometry and the calculations leading to the free data do 
not refer to the axial Killing field at all [10] (see Appendix A). In the above, we restricted 
ourselves to type II horizons primarily because they are the ones which are physically most 
interesting. On a non-extremal type HI horizon, the complex-valued seed function $a again 
captures the gauge invariant part of the free data.^ However, in absence of axi-symmetry, 
we no longer have the natural coordinate (. <I>a is now a function of two coordinates 
on A. Nonetheless, we can use the Laplacian operator q'^'^DaDb on A to decompose <I>a 
in to generalized spherical harmonics. Thus, now the geometric multipoles In,m, Ln,m are 
labelled by two integers where n labels the discrete eigenvalue of the Laplacian and m is 
the degeneracy label of eigenvectors with same eigenvalue n of the Laplacian. Should a 
physically interesting application of type HI horizons arise, it should not be difficult to 
extend this basic framework further. 

There is, however, an important subtlety. In the type II case, the presence of axi- 
symmetry led us to single out a preferred function ( and our multipoles were defined using 
the Y^{() associated with this (. As noted in section HI, these are spherical harmonics 
associated with the round 2-sphere metric q^b of (2.8). In general, they are not eigenfunc- 
tions of the Laplacian of the physical metric Qab^- Our procedure was essential to ensure that 
the angular momentum dipolc Ji agrees with the horizon angular momentum Ja defined by 
Hamiltonian methods [11]. Had we used multipoles through spherical harmonics associated 
with Qab, this equality would not have held in general. Thus the extra structure available on 
type II horizons played an important role in our definition of multipoles in sections III A and 



In the extremal case, further work is needed to understand the constraint on {qab,^a) which, in the type 
II case led to the surprising result that these fields must coincide with those on an extremal Kerr isolated 
horizon. 
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Ill C. Since it is absent on a general type III horizon, multipolcs defined using the physical 
Laplacian will not have all the physically desirable properties discussed in section III C. 

IV. DISCUSSION 

Source multipoles in Newtonian gravity characterize the way in which mass is distributed, 
higher multipoles providing a complete description of distortions, i.e., departures from 
sphericity. Multipole moments defined in this paper do the same for the horizon geome- 
try. In this sense, they can be regarded as 'source multipoles' associated with a black hole, 
distinct from the 'field multipoles' defined at infinity. In Newtonian gravity and Maxwellian 
electrodynamics, there is a simple relation between the two because both these theories are 
'Abelian': the field does not serve as its own source. In non-Abclian contexts such as Yang- 
Mills theory and general relativity, the field in the region between the source and infinity 
itself acts as an effective source. Hence one would expect there to exist two distinct sets of 
multipoles, one associated with the central gravitating body and the other associated with 
the entire system. If space-time is stationary, one would expect a simple relation between 
mass monopoles, if it is axi-symmetric, between angular momentum dipoles and if in addi- 
tion there is refiection symmetry, one would expect all odd mass multipoles and even angular 
momentum multipoles to vanish in both regimes. However, in absence of symmetries, one 
does not expect a simple relation between the two sets in the fully relativistic regimes where 
post-Newtonian and post-Minkowskian approximations fail. This point seems to have been 
glossed over in the general relativity literature on applications of multipoles. Since horizon 
multipoles were unavailable, one typically identified the field multipoles with the multipoles 
of the black hole and used them in all contexts. 

The horizon multipoles defined in this paper are likely to have three sets of applications 
which we now summarize. 

The first is to the problem of equations of motion. In the more recent literature on this 
subject, multipoles have been used in two different contexts. The first is illustrated by the 
problem of extreme mass-ratio binaries where a solar mass black hole orbits about a cen- 
tral supermassive one. For this case, there is an interesting framework [19] enabling one 
to use gravitational wave observations to 'measure' the gravitational multipole moments of 
the central object. In actual calculations, one typically assumes that the small hole is far 
from the supermassive one. One can then express the space-time metric at the position 
of the small hole as a series in inverse powers of distance from the central hole using the 
Hansen field-multipoles as coefficients, and calculate quantities of physical interest. Within 
the approximations inherent to these calculations, it is entirely appropriate to use the field 
multipoles and the framework provides a way of measuring these multipoles from the grav- 
itational wave data. However, the approximation will fail if the orbit lies in a genuinely 
strong field regime. In that case, one may be able to use the near horizon expansion of 
the metric [7]. Here expansion coefficients refer to the horizon geometry but it should be 
straightforward to recast them in terms of the horizon multipoles. 

The second context is illustrated in the study of how the gravitational field of other objects 
affects a given, small black hole (see, e.g., [20]). Calculations of the resulting distortion, 
precession, spin-flips, etc. of the black hole require the knowledge of tidal forces and torques 
exerted by other bodies on the black hole. Typically, these are calculated from interaction 
terms involving a coupling of black hole multipoles to perturbations in the Weyl curvature 
caused by other objects (at the location of the black hole). In these calculations it is 
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inappropriate to use the field multipoles of the smaU black hole since the coupling occurs 
at the location of the black hole. What one needs is multipoles associated with the horizon, 
not coefficients in the asymptotic expansion of the metric produced by the small black hole 
alone. The horizon multipoles introduced here provide the appropriate notion. 

A second set of applications is to binary black hole simulations in numerical relativity. At 
very early times when the two black holes are widely separated and at late times when the 
single black hole has settled down, the world tubes of apparent horizons are well approxi- 
mated, within numerical errors, by isolated horizons. One often wants to compare results of 
two different numerical simulations, particularly at late times. These comparisons are gen- 
erally difficult in the strong field regime because the results are tied to the coordinates used. 
Our mass and angular momentum multipoles provide a diffeomorphism invariant method 
to compare such results. Furthermore, when there are differences, they provide a tool to 
interpret them physically. 

A more interesting potential application is in the dynamical regimes, where the world 
tubes of apparent horizons can be modelled by dynamical horizons H [21]. These are nat- 
urally foliated by marginally trapped 2-spheres Sh- On each Sh, there is an obvious, well- 
defined analog of $a- Therefore, a natural strategy is to use it to define multipoles which 
would now be time- dependent. It would be very interesting to analyze the evolution of these 
moments both analytically and numerically. For example, there is a well-defined notion 
of energy flux across any Sh [21]. Is there then an analog of the quadrupole formula at 
the horizonl From the knowledge of the horizon quadrupole and its relation to the Kerr 
quadrupole, can one gain insight into the maximum amount of energy that can be emitted 
in gravitational radiation? For non-rotating black holes, early numerical simulations [22] 
introduced a notion of horizon multipoles, using tools then available. They then calculated 
a complex frequency describing the time evolution of the quadrupole, where the imaginary 
part of the frequency captures damping. They found that this frequency coincided with the 
quasi-normal frequency of the final black hole. Although the calculation of multipoles was 
somewhat imprecise, the result is simple and intriguing. Interestingly, our horizon multipoles 
provide a precise formulation of the physical ideas underlying those more intuitive notions 
of multipoles. One is therefore led to ask: Can one analyze the final stages of mergers using 
the present, diffeomorphism invariant multipoles and show in detail that the complex fre- 
quency is the same as the quasi-normal one? Such a result would provide a deep insight in 
to the relation between the highly non-linear dynamics of the horizon and the quasi-normal 
ringing associated with the linear regime. Again, the issue can be examined using analytical 
methods as well as numerical simulations. 

The final set of applications is to quantum gravity. To calculate black hole entropy from 
first principles, one needs to construct an ensemble, where the 'macroscopic parameters' 
describing the system are fixed. To be physically meaningful, these parameters have to 
be diffeomorphism invariant. For type I horizons, this is straightforward: there is only 
one gravitational parameter, which can be taken to be the horizon area (or mass). For 
type II horizons, one can not just work with mass and angular momentum because the 
horizon may be distorted by various types of hair. Even within the 4-dimensional Einstein- 
Maxwell theory where no-hair theorems hold, distortions can be caused by matter rings and 
other black holes. Even if the black hole itself is isolated, one can not automatically rely on 
uniqueness theorems and say that it must be a Kerr hole because it is physically unreasonable 
to require that the whole space-time is stationary. Since one wishes to calculate entropy of a 
black hole in equilibrium, one should only need to require that the horizon geometry is time 
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independent, not the whole universe. To incorporate physicaUy reahstic situations, then, 
one needs a diffeomorphism invariant characterization of the horizon fields. Multipoles can 
now serve as the required parameters in the construction of the ensemble. It turns out that 
they can in fact be used to calculate entropy associated with type II horizons [23]. 
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APPENDIX A: RECONSTRUCTION OF THE HORIZON GEOMETRY FROM 
THE FREE DATA 

For convenience of the reader, in this Appendix we will collect some results from [10] 
which play an important role in this paper. Specifically, we will show how the horizon 
geometry (qat, D) is determined by the free data via field equations. 

Fix a 3- manifold A, topologically x M and a 2-sphere cross-section A thereon. We wish 
to show that, given a pair of fields {qah)^a) on A, there is a non-extremal isolated horizon 
structure ([^"], qah, D), unique up to diffeomorphism, on A such that q^b is the projection of 
qab on A; Dai^ = i^d^] is the projection of Ua on A and Uai-"' 7^ 0. (We will turn to the 
extremal case at the end.) 

Introduce a function u whose level surfaces provide a foliation of A by a family of 2-spheres 
such that A is a leaf. Let Ua — —DaU and define a vector field such that i"'na^ — 1. 
Thus, is transversal to the foliation. Now, given qab on A, there is a unique tensor field 
qab on A such that : 

qabi'^0, >C,g„b = 0, and qabV'^W' ^ qabV'^W' (Al) 

for all vector fields V, W"^ tangential to A. This is the required qab- 

The D we seek is to be the derivative operator on a non-extremal isolated horizon. So, 
it must satisfy: i) Daqbc = 0, ii) [>C^, £)] = 0, and iii) Da(^ — oJa^-^ for some 1-form uja where 
I'^ijJa — is non-zero. Consider first torsion-free derivative operators D on A which are 
compatible with qab- They all have the same action on 1-forms /„ satisfying = 0, given 
by: 

Dafb=lCfqab + D[afb] (A2) 

where is any vector field on A satisfying /"'qab = ft- (While is not unique, one can check 
that the first term on the right side is insensitive to this ambiguity.) What distinguishes 
these derivative operators is only their action on Set 

Sab := DaUb (A3) 
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Then, Sab is symmetric and i'^Sab = ^''DaUb = uOa- Thus, the freedom in the definition of D 
is completely captured by the pair {ua, Sab) where Sab is the projection of Sab on the leaves 
of the foliation (defined by the projection operator Qa'' = Sa'' + riai'')- 

We now bring in the pull-backs of the field equations to A. A simple calculation yields: 

R„/^2rD[„a;6], (A4) 

where the underhne denotes the pull-back to A. Now, conditions on the stress-energy 
tensor in the definition of a non-expanding horizon imply that T^(,£" = whence R^^-^'' = 0. 
Therefore, via the Cartan identity, this part of the Einstein equation is equivalent to: 

Da{u;be')-CeUJa^O (A5) 

Now, since [Ci, D] = 0, it follows that C^uia = 0, whence Ki is constant on A. The remaining 
projection of Einstein's equations on A determines how Sab 'evolves' along A: 

A Sab = - KiSab + D{aOjb) + oja<^b " ^T^ab + ^^^a Qb'^ Red (A6) 

where D and TZab denote the derivative operator and the Ricci tensor on the leaves of the 
foliation and uja is the projection of uJa on these cross-sections. However, since [Ci, D] =0, 
the left side vanishes. Hence, if ^ 0, we have: 



Sab — — 



D{a^b) + ^J^b — l^ab + l^Qa Q.b'^ Red 



(A7) 



Thus, given the part of the stress-tensor which determines the projection Qa"^ qb^ Red of the 
space-time Ricci tensor on the leaves of the foliation, Sab is completely determined by Qab 
and oja- Since CiQab = and [C^, D] = 0, it follows that it suffices to specify {qab, oJa) only on 
one leaf of the foliation. 

To summarize, in the non-extremal case, the field equations pulled-back to A imply that 
D is completely determined by qab-, (^a, i^e and Ke is a constant on A. An examination of the 
other field equations shows that they do not constrain the horizon geometry; rather, they 
provide evolution equations off A. Thus, given {qab,i^a) on A, we can construct a triplet 
{[i],qab,D) on A = A X R providing a non-extremal isolated horizon geometry. All such 
possible triplets are diffeomorphically related because the only freedom in the construction 
is the choice of [£]. (As noted in the main text, the value of K£ is not fixed on an isolated 
horizon because of the constant rescaling freedom in [£].) 

Let us now turn to the extremal case. Analysis is identical until (A6). But, since ni 
vanishes. Sab decouples from the equation and remains completely unconstrained. However, 
we now have a constraint on the pair {qab,i^a)' 

D(ai^b) + OJaOjb - ^T^ab + \^qa qb Red = (A8) 

This equation is quite complicated and its solutions are not known in the general case. 
However, on type H horizons in the Einstein-Maxwell theory, a surprising simplification 
occurs: in this case every solution i^qabi^a) to this equation is diffeomorphically related to 
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that on the Kerr- Newman extremal horizon [16] . Hence, the freely specifiable data is just 

Sab- 
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